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We study superradiant Raman scattering from an ultra-cold, but finite temperature Bose gas in a 
harmonic trap. Numerical simulations indicate the existence of distinct timescales associated with 
the decoherence of the condensed versus thermal fractions, and the concomitant preferred scattering 
from atoms in low lying trap states in the regime where superradiance takes place on a timescale 
comparable to an inverse trap frequency. As a consequence the scattered atoms experience a modest 
reduction in temperature as compared to the unscattered atoms. 



I. INTRODUCTION 

Bosc-Einstcin condensation (BEC) is characterized 
by the appearance ofF-diagonal long-range order, with 
phase correlations between points spatially separated 
over macroscopic distances. These correlations are 
quantified by the first-order spatial correlation function 
G(i)(x,x') = (^/'t(x)v;(x')), where ^^■^(x) is the field oper- 
ator that creates a particle at position x. In the thermo- 
dynamic limit G^^^(x, x') becomes finite-valued for in- 
finitely separated spatial points below the critical tem- 
perature. In ultracold atomic systems the presence of 
long-range phase correlations is experimentally manifest 
through the direct observation of matter-wave interfer- 
ence patterns. Hence, matter-wave interference has been 
used as a diagnostic tool to demonstrate the appear- 
ance of the superfluid to Mott-insulator transition [l|, 
the growth of spatial correlations during the formation 
of a Bose-condensate from a non-equilibrium situation 
after sudden quenching across the transition point 0, 
and to measure the critical exponent characterizing the 
divergence of the correlation length near that point Q . 

Another tool recently employed to probe the BEC 
transition is the selective sensitivity of superradiant scat- 
tering to the condensed fraction of an ultra-cold Bose 
gas beneath the critical point for condensation Q . Early 
demonstrations of superradiance were realized in thermal 
gases at temperatures in the several hundreds of degrees 
range, @, @, 0, Si- The majority of these experiments in- 
volving initially inverted atomic transitions in relatively 
high-density sanrples, a situation often referred to as su- 
perfiuorescence [H, Q. By contrast, recent experiments in 
ultra-cold atomic vapors involved off-resonant light scat- 
tering in very low density samples. This regime was first 
studied in 1999 and has since been the focus of several 
experiments, such as in the context of coherent matter- 
wave amplification [ll|, [l^, The majority of these 
studies used almost pure condensates in which superra- 
diant enhancement is strong due to the near absence of 
Doppler dephasing. In addition to related theoretical [l3| 
and experimental work on coherent atomic recoil lasing 



0, [3| , to our knowledge only two recent experimental 
studies considered in detail superradiance from ultracold 
thermal vapors and its behavior as a function of temper- 
ature 

One important feature of superradiance from a bosonic 
atomic vapor cooled below the BEC critical tempera- 
ture is the existence of two well-separated decay time 
scales, the fast one attributed to the Doppler dephasing 
of the thermal fraction of the gas, and the other, slower 
by roughly an order of magnitude, due to the condensed 
fraction. These separate time scales were observed by 
Yoshikawa and co-workers [l^ and exploited by Sadler 
et al. Q to selectively image the condensed fraction of 
the vapor over a range of temperatures above and below 
the critical temperature. In the latter experiment [3| su- 
perradiance enhanced absorption of light scattering from 
a prolate sample of condensed bosonic atoms was mon- 
itored in real time during the superradiance process. A 
main goal of the present paper is to provide a theoretical 
description of superradiance in ultracold bosonic gases 
to explain the key temperature dependent aspects of the 
above-mentioned experiments [1, [l3l , that enable its use 
as a probe of coherence in Bose condensates. As such, 
the present paper is an extension to finite temperatures 
of our previous study [3 of the spatial inhomogeneity 
observed in absorption images of the probe light at T = 0. 

The paper is organized as follows. Section II discusses 
our theoretical model of a non-interacting, ultra-cold 
atomic gas undergoing Raman transitions in a spheri- 
cal harmonic trap. Section III presents numerical results 
demonstrating the existence of two decay times in the su- 
perradiant signal. Section IV shows that in that regime 
where the superradiance time scale becomes comparable 
to a trap frequency, superradiant scattering occurs pre- 
dominantly from atoms in the low-lying energy states of 
the trap, hence it may be used to selectively probe those 
states. We also remark that superradiance may leave the 
'cold' and 'hot' atoms in two different electronic states. 
Finally, section V is a summary and conclusion. 



II. FINITE TEMPERATURE THEORY 
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of Standards and Technology, Boulder, CO, 80305. 



We consider TV non-interacting ultracold atoms in a 
spherical harmonic trap. Fig. 1(a). The atoms undergo 
A-type Raman scattering between two electronic ground 



2 




FIG. 1: (a) Experimental setup - a spherically shaped atomic vapor is driven by a laser beam incident along the s-axis. The 
end-fire modes travelling along the z-axis are superradiantly amplified. The pump light is assumed to drive a A-type Raman 
transition for which uol ^ i02 ^ 5l, see Fig. 1 (b) 



states |1) and |2) via an excited state |e), see Fig. 1(b). 
We assume that the transition |1) — > |e) is driven by an 
off- resonant classical pump laser Ei(i) of frequency oj^ 
propagating along the x-axis of the trap. For a spherical 
trap there are no geometrical effects [I3 that lead to the 
selection of particular superradiant modes as would be 
the case for a cigar-shaped trap, but a preferential direc- 
tion can be selected by placing the sample in a low finesse 
ring cavity or by first applying a Bragg pulse to generate 
a seed matter wave [111, [l^l . We assume that this is the 
case here and that as a result the dominant superradi- 
ant modes propagate along the z-axis. We refer to them 
as end-fire modes in analogy to the case of elongated 
samples. The polarization of the pump light is also cho- 
sen parallel to the z-axis. For that polarization Rayleigh 
scattering into the end-fire modes is suppressed due to the 
angular radiation pattern of Rayleigh scattering. It does 
however allow Raman scattering into circularly polarized 
end-fire modes which we describe quantum-mechanically 
in terms of bosonic creation and annihilation operators 
djj. and ak, respectively. The total electric field is then 



E 



TlUJi. \ 



2enV 



akWe*" + /i.c. 



where the incident laser field envelope is taken as 
constant in amplitude. In terms of the detuning 



we have 



■ IjJl 



5l. 



(1) 



(2) 



We proceed by introducing bosonic matter-field cre- 
ation and annihilation operators ^\ (x, t) and ipi (x, t) , 
that create and annihilate, respectively, an atom at time 
t and position x in electronic state \i) = |1), |e) or |2), 
with 



7/.,(x,t),V;Vx',i) =%<5(X-X'). 



(3) 



Taking wi = the Hamiltonian of the atom-field system 
is H = Ho + He, with 



Ho = fji^kaj 



E 

i=l,2,e 



!x |?ia;e'0e(x)'0e(x) + ^W2V'2 WV'2(x)| 



2m 



^.(x),(4) 



while 

Hr = 



dx{E-d V)t(x)7/'i(x) + V'|(x)V;2(x) -t-/i.c.} 



(5) 

describes the electric dipole interaction between the 
atoms and the electromagnetic field, d being the dipole 
moment which we take to have the same magnitude for 
both transitions. 

For large enough detunings the excited state is not sig- 
nificantly occupied and may be adiabatically eliminated. 
Introducing slowly varying operators ipi = ^^e*'^'* and 
Ok = Q.ke*'^''* and using the rotating wave approximation 
yields then the effective interaction Hamiltonian 
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(6) 



Here fti = dEL/h is the Rabi frequency of the incident 
field and flk — d^Jukl {2eQhV). 

Expanding the atomic field operators in states of the 
trap, 



(7) 



1=0 



where iJo</'n(x) = hi>n4'n{'^) and the subscript n labels 
generically excitations in all three dimensions of the trap, 
n = {^x, ny,nz}, and = {nx+riy +71^ + 3 /2)ujt, substi- 
tuting the expansion ([7|) into Hamiltonians ([6]) and ([4]), 
and denoting the expectation values (at) ~ Ok, wc obtain 
the Heiscnberg equations of motion 



dr 



dr 
dr 



k 



km 



(10) 



r 



Here the dimensionless time t = flut where = 
^iL^k/^L is an effective 'two-photon' Rabi frequency, 
j/j'j ~ Vn/O.R, the slowly varying operators Cni — Cnie^^^''' 
and the overlap function 



'7nj(q) 



,(x)(/.,(x)e^'i-dx. 



(11) 



where the recoil momentum is q = — k. A grayscale 
rendering of the absolute value of rjnj (q) is plotted for the 
first 25 trap levels and the trap and light field parameters 
introduced in section HI. Finally, the time evolution of 
the end-fire modes is given by 

9ak 



ij] W (q) (42Sni)e^('''"-^")* - rak, (12) 



where the last term is a phenomenological decay added 
to account for the escape of the photons from the sam- 
ple. This term should in principle be accompanied by 
appropriate quantum noise operators to guarantee that 
the commutation relations of the operators Ok are pre- 
served at all times. These noise operators are also impor- 
tant in that they trigger the superradiant amplification 
and determine the field fluctuations in its early stages 
[13, HH- This paper concentrates however on the later, 
classical stages of the process, and we follow the standard 



procedure of introducing a small classical seed for the 
initial amplitudes of the cnd-firc modes, which are then 
described as classical fields. At that level of approxima- 
tion the phenomenological decay term in Eq. (jl2p is ap- 
propriate. Note also that by factorizing the expectation 
values of the optical and matter-wave fields, we have ne- 
glected any quantum correlation and entanglement that 
may build up between the optical and matter-wave fields 
fields. 

The atoms obey initially an equilibrium Bose distribu- 
tion. Wc determine the distribution for excitations in the 
z-direction by averaging over the x and y-directions, 



(cllCnl) 



1 



ahu:t{nx+ny+nz)/kBT ^ ' 



(13) 



where A/" is a normalization factor chosen to fix the ther- 

oo 

mal fraction of atoms, rith ~ (siiCni), at the value 

rith = N{T /TcY determined by the temperature of the 
sample [22|, where Tc is the critical temperature. 

We conclude this section by remarking that in order 
to keep the computations manageable we have assumed 
that the light field intensity is uniform within the sam- 
ple, thereby neglecting spatial effects that have previ- 
ously been shown to play an important role in the growth 
of the cnd-firc modes [II, HI, [2J . The role of these effects 
at finite temperature remains therefore an open question. 



FIG. 2: Gray scale rendering of the absolute value of the 
overlap function |7?mn(fc)| for the first 25 trap levels. Here 
we've assumed a trap frequency ut = 2n ■ 800 rad/s and fc_L ~ 
27r/795 nm corresponding to the _F = 1 — > _F' = 1 transition 



In the following simulations we also assume for simplic- 
ity that the atoms recoil only along the z-axis and we are 
therefore only interested in excitations along that direc- 
tion. 



III. DEPHASING 

This section summarizes numerical results that illus- 
trate the dynamics of the thermal and condensed frac- 
tions of the vapor, showing that these dynamics are char- 
acterized by two distinct time scales. Specifically, we sim- 
ulate a pump-probe technique that was experimentally 
implemented by Yoshikawa and co-workers [TtI to mea- 
sure superradiant coherence. This procedure entails ap- 
plying a superradiant pump pulse until the end-fire mode 
intensity reaches a maximum. The growth of the end-fire 
mode occurs in response to the build-up of atomic polar- 
ization (a polarization grating forms), Eq. (jl2|l . At that 
point the pulse is turned off for a variable delay time, t^, 
during which the polarization grating deteriorates due to 
Doppler dephasing. The pump beam is then turned back 
on to its original intensity to probe the remaining coher- 
ence. Had the polarization grating remained intact the 
superradiant intensity would immediately return to its 
pre-delay level, however, due to Doppler dephasing the 
post-delay intensity is reduced. The ratio of the pre- and 
post-delay end-fire mode intensities, 

C = |ak(Tmax + Td)|V|ak( )r, (14) 

is therefore a measure of the decoherence of the superra- 
diant pulse. 

Figure [3)Ja) illustrates the detector response (upper 
curve) to a square pulse pump-probe sequence (lower 




FIG. 3: (a) Pump-probe spectroscopy. The lower curve shows 
the pump-probe sequence and the upper curve the resulting 
end-fire mode intensity. The ratio of the peak heights di- 
rectly after to directly before the delay is used as a measure 
of coherence, (b) Dephasing. Dotted line: post- to pre-delay 
intensity signal ratio for a vapor with 0.95 condensed fraction. 
For a condensed fraction of only 0.1 (solid line) the decay is 
bimodal, the initial rapid decay being attributed to the ther- 
mal fraction. Dashed line: intermediate case, ric = 0.5. (In- 
set) Solid circles: polarization | (cJ„2'^oi) | ; open circles Fourier 
components of the signal for — 0.95 in the main figure. 
The dotted line envelope shows the overlap function l?7mo(q)|. 



curve). Here the end- fire mode intensity was convoluted 
with the response function of the detector, assumed to 
have a (dimcnsionlcss) response time of r' ~ 0.01. 

Our simulations are for a sample of = 10^ atoms in 
a three-dimensional, isotropic harmonic trap with trap 
frequency = 27r • 800 rad/s and, with reference to [j| 
we use the = 27r/795 nm corresponding to the F = 
1 F' = 1 transition of ^^'^Rb. We use ujt/flR — 0.5 and 
included 400 trap levels in the simulations. A tight trap 
was chosen for computational reasons: For that value of 
uif the occupation of levels higher than n = 400 remains 
neg ligible at all times for the conditions that we consider 
[2^. We choose ri^/Jlfc = 1000 and neglect laser pump 
depletion. Since the decay of the light field due to its 
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FIG. 4: Summary of numerical simulations for ric = 0.5. (a) Circles: Ratio of the total number of thermal atoms (n^^ ) to 
total number of condensed atoms {ric ) transferred from the state |1) at peak superradiant intensity ; triangles: fraction of 
total number of atoms scattered to the state |2) (condensed and thermal included), as a function of ujt/^R- (b) Fraction of 
atoms left in ground state |1) at the end of the pulse complete by trap level for LOt/O.R = 0.25, 0.5, 0.75 and 1.0. (c) Initial 
distribution of atoms by trap level (solid circles) and distribution of atoms in state |2) at the end of the pulse for uot/^R ~ 1-0. 
The initial distribution was multiplied by a factor of 10 for visibility. 



escape from the atomic sample, see Eq. (fT2|l . is the fastest 
process by several orders of magnitude the state of the 
scattered light field is nearly instantaneously determined 
by the state of the atomic fields, so that 

flk « ^ E'?'"" (42C«i)e^(''"-''")*. (15) 

mn 

Either a seed atomic polarization or a seed end-fire 
mode occupation is required to initiate the superradiant 
growth. In this paper we usually assume that the end-fire 
mode has an initial value Ok = a/TO and that all atoms 
are in the state |1) obeying the distribution (fT3|l . The 
first numerical iteration of Eqs. ([5|)- (|10|) then creates a 
small initial polarization. Thereafter we use Eq. (jlSp to 
determine the instantaneous value of the end-fire mode, 
with r w (10^ ~ 10^) X riR. 

Figure |3)Jb) plots ^{rd), the post- to pre-delay inten- 
sity ratio, Eq. (|14[) . as a function of the pump-probe de- 
lay time for three different sample temperatures. Since 
the supcrradiancc process is initiated by quantum fluc- 
tuations leading to spontaneous emission of photons into 
the end-fire modes, shot-to-shot fluctuations occur in the 
relative amplitudes and phases of the left- vs. right- 
propagating end- fire modes [l^. To account for these 
fluctuations we chose, for the purposes of Fig. [3]Jb), the 
initial seed amplitudes from a random Gaussian distribu- 
tion centered around at = \/To with a standard deviation 
of cr = \/T0/2 and with a random overall phase. Each 
curve in Fig. represents an average over fifty runs of 
pump-probe experiments. 

The dotted line corresponds to a nearly pure conden- 
sate, with a condensed fraction iic = 0.95. After an initial 
increase for short times due to the continued build-up of 
the atomic polarization , the signal ratio decreases to half 
the initial intensity with a decay time of Tc ~ 0.3. 

The slow oscillations in ^{rd) can be understood by 
comparing the Fourier transform in one single shot of the 
experiment to the components of the atomic polarization 
l('^m2'^oi)| between the trap ground state and its excited 



levels, see insert of Fig.[3l^b). As expected, there is a one- 
to-one correspondence between these frequencies. Hence 
the oscillations in ^(r^) are a signature of the specific trap 
characteristics. In the inset the atomic coherence is only 
significant between even trap states. This is a combined 
consequence of the two end-fire mode amplitudes and 
phases having been chosen equal in this shot, and of the 
form of the overlap function ('7jo('?)) between the lowest 
trap level and higher levels, which are alternately purely 
real and purely imaginary, due to the harmonic oscillator 
states being alternately even and odd. In Eq. PU)) . only 
the last term on the right-hand side contributes initially. 
Then, by choosing the initial overall phase of the left- 
propagating end-fire mode to be the same as that of the 
right-propagating one, the left and right contributions for 
all odd trap levels cancel out since the overlap functions 
for those modes are purely imaginary. The cancellation is 
no longer exact if either the phases of the left and right- 
propagating end-fire modes are not the same or if their 
initial amplitudes are not equal, as is in general the case 
due to fluctuations. 



As the temperature approaches Tc, £,{Td) undergoes a 
rapid initial decay on a timescale rth roughly an order of 
magnitude faster than the slow decay characterized by 
Tc- This is illustrated for Uc = 0.1 as the solid line in 
Fig. [Sjb). These two time scales were experimentally ob- 
served by Yoshikawa et at, see Fig. 3 in Ref. [l3|. The 
appearance of a slow decay below the critical tempera- 
ture Tc, combined with the disappearance of the rapid 
decay at T — > 0, clearly points to the fact that they are 
associated with the condensed and thermal fractions re- 
spectively. The existence of these two time scales was 
exploited by Sadler et al. [1| to image selectively the 
condensed fraction of an ultra-cold Bose gas below the 
critical temperature while remaining blind to the ther- 
mal fraction, as we discuss next. 
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IV. THERMALLY SELECTIVE SCATTERING 

This section further discusses how the Dopplcr broad- 
ening of the thermal fraction inhibits superradiance when 
the effective Rabi frequency fi^ becomes comparable to 
the trap frequency, thereby making superradiance a sen- 
sitive probe of the condensed fraction of the sample Q . 

The circles in Fig. IHa) show the ratio /ric , the 
number of atoms transferred to the electronic state |2) 
from the excited trap levels to the number of atoms trans- 
ferred from the trap ground state to |2), at peak end- fire 
mode intensity and as a function of ujt/^R for an initial 
condensed fraction ric = 1/2. As the ratio uJt/^R in- 
creases, fewer thermal atoms participate in the superra- 
diant scattering as a consequence of Dopplcr dephasing. 
The triangles, which show the total number of scattered 
atoms, confirm that indeed, the bulk of the superradiant 
emission is associated to the condensed atoms in that 
case. 

This behavior is further illustrated by plotting the frac- 
tion of atoms left in the electronic state |1) at the end 
of the superradiant emission as a function of trap level. 
This is illustrated in Fig. [4{b) for the cases (from lowest 
to highest curve) uot/^R = 0.25, 0.5 and 0.75 and 1.0. 
When superradiance occurs on a time scale faster than 
the inverse trap frequency, all levels contribute roughly in 
proportion to their initial occupation [26j . But as uit/^R 
becomes of order unity, the lower trap levels contribute 
significantly more in proportion to their initial occupa- 
tion than the higher levels which Doppler dephase faster. 

The selective scattering of atoms initially in the deep- 
est trap levels implies that the scattered atoms are in 
effect colder than the unscattered atoms — a situation 
somewhat akin to evaporative cooling, except that in the 
present case it is the warmer atoms that are "left behind." 
The squares in Fig. \M,c) show the atomic distribution of 
atoms in electronic state |2) at the end of the superra- 
diant emission for the first 20 trap levels for ric = 0.5, 
ujt/flR = 1.0 and N = 10^. The solid circles show the 
initial distribution of thermal atoms in state |1) (the oc- 
cupation number of the lowest trap level is off scale and is 
not shown here and the initial distribution was multiplied 
by a factor of 10 for visibility). 

We use the root-mean-square deviation of the center- 
of-mass energy of the atoms as a measure of their temper- 
ature, and compare its initial value Aei for atoms in the 
state |1), to that of the atoms in the electronic ground 
state 1 2), Ae2, at the end of the superradiant emission. 
For LOt/VtR = 1.0 we find Ae2 = Aei/2, a modest reduc- 
tion in rms energy. This is despite the fact that several 
trap levels wind up being significantly populated, as op- 
posed to only the lowest trap level being significantly 
populated initially and follows from the comparatively 
few thermal atoms undergoing Raman scattering, as we 
have seen. For the case at hand, 58% of the particles 
in the state |1) initially occupy the first 10 trap levels, 
but after scattering 82% of the particles scattered to |2) 
occupy these 10 levels. 



This effect can be enhanced if the trap frequency is 
chosen so as to approximately match the photon recoil 
frequency, Tiq^ /2M ^ tot, the spatial width of the ground 
state of the atom trap thus being of the same order as the 
wavelength of light. This im plies that the Lamb-Dicke 
parameter, 77 = kL\/'h/2mujt, which serves as a measure 
of the amount of coupling between the motional and in- 
ternal states of the atoms, is of order ~ 1. The overlap 
function associated with the ground state, r?om(q), is then 
appreciable only for m = 1, and the condensed atoms will 
only recoil to that level, instead of several levels as was 
the case in Fig. [DJc) . 



V. CONCLUSION 



In this paper we have explored theoretically various 
aspects of superradiance from ultracold, but finite tem- 
perature atomic gases. In particular we've illustrated 
explicitly the existence of two well separated coherence 
timescales associated with the condensed and thermal 
fractions of the gas respectively. In the regime where 
Doppler dephasing plays an important role, we found 
that superradiant scattering takes place predominantly 
from atoms in the lowest lying trap levels as these are 
less prone to the Doppler dephasing. That effect was ex- 
ploited in [3] to image selectively the condensed fraction 
of an ultra-cold gas. We demonstrated that as a conse- 
quence atoms in the scattered state have a reduced rms 
deviation from the average trap energy as compared to 
the atoms in the initial state, and are in that sense colder 
than atoms in the initial state. 

Our calculations considered a system with end-fire 
modes recoiling perpendicularly to the pump beam. In 
a system where the pump beam propagates collinearly 
to the scattered beams, such as an elongated condensate 
pumped along the long axis, it is possible to suppress 
Raman scattering into the backward end- fire mode, as 
compared to the forward end-fire mode, by exploiting 
the same dephasing effects that lead to thermally selec- 
tive scattering. This results because the recoil momen- 
tum associated with the forward end-fire mode is small 
compared to that due to backward scattering. Hence it is 
possible to realize conditions such that the forward scat- 
tering Lamb-Dicke parameter rjf <^ I and the backward 
scattering one rji, ^ 1, implying negligible coupling be- 
tween the ground state and higher motional states for 
forward scattering, but strong coupling and accompany- 
ing dephasing for the backward scattering case. This 
effect would not be present in Rayleigh scattering exper- 
iments using the same geometry, since there the pump 
beam simply sees a phase shift in the forward direction 
and only the backward mode is superradiantly enhanced. 
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